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CHAPTER 1 


INTRODUCTION 


The statistical model developed and stated in NAVWEPS 
OD 29304 "Guide Manual For Reliability Measurement Program" 
was developed for the Department of the Navy, Special 
Projects Office. It contains very versatile procedures for 
obtaining both system reliability point estimates and system 
reliability lower confidence limits using a variety of test 
data on component parts, subassemblies, assemblies, and 
subsystems. 

This document makes the following assumptions on the 


components of the series system; 


(1) Constant failure rate 
(2) Additivity of stress effects 
(3) Independence of part failure 


(4) Failure rate constancy 


The purpose of this thesis is to check the accuracy of 
this statistical model without the assumption of a constant 
failure rate. The model will be checked for a system conm- 
posed of four components in series, using various sample 
sizes, mission times, and combinations of the gamma and 
truncated normal failure distributions, The statistical 
model will be analyzed by computer simulation using a Monte 
Carlo method to generate the gamma and truncated normal test 


times. The parameters of the distributions were selected so 


that the failure rate curves would be representative of 
those curves used by contractors or that the System relia- 


bility was a preselected value. 
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CHAPTER II 
THE STATISTICAL MODEL 


Section three of reference 1 cContaims a detailed 
explanation of the statistical model which is being analyzed. 
The following definitions are applicable parts of the model 
used in the computer simulation for the lower confidence 
limit of a series system with four components. For a given 


Sample size, N 1=1,2,3,4 and a predetermined test time, 


i? 
Toi? i1=1,2,3,4 a best estimate of the failure rate for 


component i is given by: 
i_ ed (1) 


N. 

i 
Where S4 = weyo2d? 
accumulated on the N components of type i, and fs is the 


is the sum of all test times 


number of failure times which are less than TS From these 


ee 
best estimates of failure rates, the system reliability 


estimates are given by: 


> 
It 
0) 


(2) 


éN K 
where ho = ye (3) 


The model then uses the following formulae in order to 


eompute the bower confidence limit onesystem reliability: 


did 


| A * fins 
oe + (Hepes + \ Wie (Ka)°C ++ (xa) te 





‘A 
OO (4a) 
k 
ists we f, > @ 
i=1 
ac 
* _ tks’ = (4b) 
u n joo ae 
k 
if £ f, = 0 
i=l 


where n is the number of tests. 
Here Net the upper confidence limit of the failure 
rate for a system with k units in series, and K is a per- 


centage point of the normal distribution. § is a correction 


wk X 


factor for small values of ay and C = oe All times 


> 
iH 
-~ 

=| 
be 


are in mission units. 


A 
Reference 1 then defines R to be tne Lower con- 


s,L(a) 
fidence limit on system reliability; 


A 
A -A, 
As, L(a) Sa ae (55) 
It must be noted that since equation (2) is only a best 
estimate for the failure rate, the statistical model ex- 
plained serves only as an approximation for system relia- 


loa dbabe see 


ip 


CHAP CR Lit 
SIMULATION PROCEDURE 


The system to be simulated will consist of four com- 


ponents in a logical series, with a system reliability, Ro» 


which may be expressed as: 


wnere R, = P(T>1) is the true reliability of the i°) com- 


ponent at mission time one. It is desired to find a lower 


100{(1-a)% confidence limit for R., namely the statistic 


A 


R Be El is in fact the exact 100(1-a)% lower 


s,L(a)’ Sialic) 
GCenridéence Limit for Ro» then: 


=o (6) 


holds. 
From equation (6), R, should be the (1-a) percentile 


point of the simulated distribution of R L(a)° 
9 





For the investigation a will be chosen as .20 and Ay a.) 
Will We deGined as they (1 —a)epencentidespoint sof sbhe- diieitxi— 


“NN 
bution of R. L(a)° In order to investigate the distribution 
9 


-= 


of the random variable, a simulation by digital computer is 
used to generate the values needed. Five hundred values of 


A 
ie L(a) will be calculated and ordered such that: 
9 


aN A A IN 
eee eyes Licyee # Mere y600) 


Sauce Ghemvelue=or @ chosen for The simulation is..20,.09e 


ox 
value of R will be found such that 20% of the ordered 


s,L(a) 

values are above it. Thus the 80°" percentile of the distri- 
a 

Duron Of As sL(a)? Ada)? is found and should equal Ro. 

The measure of accuracy will therefore be the value 


[R, - The accuracy of the method will be checked 


for different combinations of the failure rate distributions 
among the four components. The sample size, Ny» to be tested 
Will consist of 50, 10, and 500 units. The mission times 


for each component, Toa? will bes: 


(1) .5 mission units 


(2) 5.0 mission units 


(3) T, such that ~ ge tala] <j (7) 
£(1) 
where 2(1) = 3774 (8) 


The third test was selected as an ad hoc rule because of the 
intuitive feeling that tnis time would give a good measure 

NN 
of accuracy to the procedure. The distribution of Rs tla) 


Will therefore be dependent upon: 
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(1) a 


(2) Ny 1 = dhe 2 oe 

(3) Fg 1 ee 2 ee 
Ol 

(4) A, 1 = is Cie sit 


The mean and the standard deviation of the distribution 
ran 
of Rs ,L(a) Will also be calculated in order to better eval- 
uate the accuracy of the model. The procedure for the 


computer Simulation of the model is: 


(1) Generate random times, Ty, 1 = 6 2eee f= en 


for each of the four components in series. 


| ¢ t t 
Tay Typ M5 2 2 8 May 
t 
‘ ’ ' ‘ 
a al le 32 ah 33 i. : we Ih 3N 


aaa T" yo sett e e ° Pay 


(2) Test these times against the selected mission time, 


uh i=1,2,3,4 and if T* » record a failure 


< 
One. 5 al Toa 
or if T a3 > To,» record a mission success and then 


set 


oe Cia: . a 7 a 
i ij ge es To4 | 


id Meg iy. - | 
Ort ehl 


(3) Brommthenmequat donee 1.) e—miG>)qeu lige! Aaechapter II, 


i 


{\ 
calculate the resulting values of eke L(a} and Ho 
9 


(4) Repeat the above steps five hundred times. 


ils 


(5) Order the computed values of ne 
(6) Calculate the mean and standard deviation of the 
resulting distribution and find the 80%) per- 
Zemwanpe point, A(, -a)° 
The complete computer simulation is found in Appendix C. The 


procedure used is shown above and comment statements have 


been used to aid in following the steps of the progran. 
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CHAPTER IV 


RESULTS 


For the procedure stated in Chapter III, the failure 
times generated will be either gamma distributed or trun- 
cated normal distributed. The parameters of the distri- 
butions were chosen so that the reliability of the component 
Was a preselected value or the failure rate curves would be 


representative of those curves used by contractors. The 


following cases were simulated: 


CASE COMPONENT DISTRIBUTION PARAMETERS ie, 
il We gamma r=1.5, s=.036 .9950 

3,4 tr. normal A=11.28, T=4 .9950 

1a ie ganna r=1.5, s=.036 9950 

3,4 CL, norme lL AM=26, ¢ =L0 9984 

iat Lae gamma r=1.5, s=.036 ©9950 
3,4 tr. normal w=-30, C=45 9707 

Tay eRe Cheat gamma r=1.5, s=.036 .9950 

V eo ay tr. normal A=8, J =40 9810 

VI a oH tr. normal AA=20, FT =40 ~9850 

Vil 1 ey tr. normal AA==—30, THH5 ~9707 


Bach case is specified by the type of distribution and 
parameters applied to the four components in series. The 
results in Table I are for the first three cases using all 
possible combinations of sample sizes and test times. fThe 
mesults of cases IV, V, Vl, Sma VII are @iven in Table II. 
These cases are different from the first three in that mis- 
cellaneous test times and sample sizes are used in the simu- 
lation, 


Ie 


The first five columns in the tables are self explana- 
Cory and tne last is the value of TT which is represented by 
ra) 
Gao smongiyoulbes. IMe eee! hs Toy- It is a measure of tne 
amouacl Of Seesting, nequimed, reskatiwe to the component un- 
1ef=) hoM=Venk ls olieats- 
Appendix A contains discussions of the gamma and 


truncated normal failure distributions. Appendix B contains 


the computation of the ad hoe test times used in Table J. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 


A partial explanation of the procedure developed in 
reference 1., which is used to compute an upper limit on 
failure rate estimates and lower limits on reliability 
estimates, has been presented in Chapter II. A simulation 
of this procedure, without the assumption of an exponential 
failure rate, was proposed in Chapter III and the results 
tabulated in Tables I and II. 

Examination of the results of the first three cases 
shows that the average accuracy, D, for the various planned 


test times, PTT, are as follows: 





Many more simulations would have to be conducted, using 
Various parameters of thé failure rate distributions, in 
order to establish final results. But from the cases 
presented it would seem that the test time selected from the 
ad hoc rule stated in Chapter III, produced the best accuracy 
regardless of sample size. From Cases IV and VII the average 
accuracy using the ad hoc rule for test times was .0085. In 
Cases V and VI téSt tim@e of 1.0, 5.0, 10.0, 15.0, and 20me 
mission units were used and the accuracy varied from a low of 


~0000 to a high of .0071. 
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APPENDIX A 


A. Gamma Distribution 





iguanas oan - St 

f(t) = Fate Ti (9) 
r-1 - st 

Z(t) = sts = (10) 
e “st , (a) 
sO es 


To evaluate the parameters and meet the requirements that 


Re = .995 the following procedure Was used; 
ie) erat | 
a, = Paek.) 
oe PLT >1 |] 
R, = p[2sT > 2s | 
| Geer 
By - PIAS. > 2s] 
995 = PIXS.. > 2s] 


Therefore for a given value of tne parameter r, and the con- 
ponent reliability, a chi-square table look up would give the 
value of the second parameter s. For r= 1.5, s = .036. The 
specific value of r was chosen because the failure rate curve 
is ee of those curves used by contractors. A 
plot of equation (10) with the selected parameters is shown 


in Figures IA through IIIA. 
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B. Truncated Normal Disceribution 




















+ (ian 
b o sa 
(CE = — = ale (u ) (11) 
2) ove ? 
where u= vaca 
— 
1 
- = EZ 
: @ () a 
Lite 
GC) = ———e 
Z(t) RCE 
a4 (EA 
- bg = 
a(t) = Ove (13) 
b ak (ae 
64 Gis dt 
C720 e 
iL 
ata) - Pw) 
1- @(%) 


mer Case I, BR, Was specified at .995 and the following pro- 


cedure was used to evaluate the parameters GQand “™. 


hee Pe Le 
i 


= ‘:—p LM La 
= = Aaj 
A, =) ei (il o( Hh) 


- a (_ LAA 
£995 = v(Q ( =a)) 
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In Case I, O and AAwere selected to be 4 and 11.28. From 
equation (12), b= 1.0024. For Cases II and III the values 
of the parameters were chosen such that tne failure rate 
Cummes WOUldwpe Llatter and again better approximate the 
failure rate curves used by contractors. In Case III tne 
restriction of a lower component reliability was also placed 
on the parameters. Parameters for Cases V and VI were 
selected as miscellaneous values in order to give more con- 
plete results to the analysis. Plots of equation (13) with 


the selected parameters are shown in Figure IA through IVA. 
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APPENDIX B 
COMPUTATION OF AD HOC TEST TIMES 


The test times for Cases I, II, and III were selected as 


follows: 


wy 5 misstonsunits 
2s 5. Missionsunits 


ee: al 
3. Ty, such that iia Z(1) (7) 


The computation of the planned test times under the ad hoc 


rulestfor the distributions used isegiven below, 
A. Gamma Distribution 
R(t = PC tb) 


substituting into equation (7) 


ilk , “ls 
—~—— In(R(T..)j=- —— ln(P(T>T, )) 
Ts OL - oi 
ab 
= —- —— In(P(2eT>2sT_}) 
T. oL 
One 
i eee é 
- tT in(R(T,, )) “= eam nen In(P(X5,.>2sT., )) 


Therefore the following equation must be solved for Toa? 


1 2 
= i In(P(XS > 2sT. 


)) = 2(2) 
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Solving equation (10) with the parameters r = 1.5 and 


s= .036, Z(1)=.0076. Substituting anbouticwmovesc uae: 


- in (P(X5 > .072 T,,)) = .0076 T,, 


Solving gives T., = 2.4 mission units. 


i 
B. Truncated Normal Distribution 


ett) eat leet) 


substituting into equation (8): 


a -— 4 
= ao in(R(T,, )) =- 7 in (P(T>T,,)) 
Ol et 
- 4 in (R(t,,)) = - Le in(o(1 -$ (PQ4=A)) 
Toi ial Oi wT 


Therefore the following equation must be solved for T6542 


T . = 
- in (o(1-$(-224))) = 202) Ty (144) 


Solving equation (13) with the parameters for Case I, 
Z(1) = .0037. Substituting into the above equation: 


eS 


iL 
- ln (1.0024 (1 BF kaa 


i oe . 00s aa 


maeretore for the truncated normal distriburion) in Case i, 


T .=1.8 mission units. Solving equation (14) with the 
oi 


515. 


parameters for Cases II and III, the following times result: 


Case ITI a” =liwe> mission units 


‘ 





een MODIS (C 


The gamma and truncated normal failure times were com- 
puted using a Monte Carlo method. A pseudo uniform random 
number generator was used to generate a normal (0,1) random 
variate, V. The random number generator and the compiling 
of the normal variate was programmed in machine language 
and was used as a subroutine function called SNORM. The 
random number generator used is called URN and is available 
at the N.P.G.S. computer facility library. 

The gamma distributed failure time is compiled from a 
cniesquare value with er degrees of freedom and a parameter, 
S, such that the gamma variate = a Sinee ea is equal 


to the sum of 2r squared normal (0,1) variates, the gamma 


distributed time is compiled in the computer simulation by 


2r 
2 V 


Boe Value, ae The normal distributed failure times 


2 


are truncated at zero and they are computed from the gener- 
ated normal variates. They are set equal to QO V+A, 
where 6 and y are the parameters of the distribution. 

The general flow of the simulation is to test each 
Component seperately for N, replications, The eCstimave or 
that component failure rate, \,, is given by statement 
number 586, which represents equation (1). After each com- 
ponent has been tested N, times and it's failure rate 


An 
estimated, the system failure rate is given by Ag. Tne 


pommy CStimave Of sysvem reliability is then compuvced an 
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statement 591, which represents €yuation (2). The number 

CRs tilmosmanatLne standard deviations for the trials of 

all four components are compiled in statements 587, and 588. 
ifemrecsurcrme trOrmuliac for the lower confidence amits of 

the failure rates are given in statements 616 and 624 de- 
pending on the number of failures, as explained in equations 
(Ya) and (4b). The whole process is repeated five hundred 
times and the distribution ordered. The mean and standard 
fevLaciton 10or Rea and R. the reliability point estimate, 


are also calculated in tne computer simulation. 
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